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Based on a model Hamiltonian with a d-wave pairing interaction and a competing antiferromag- 
netic interaction, we numerically study the site dependence of the nuclear spin resonance (NMR) 
relaxation rate T^^ as a function of temperature for a d-wave superconductor (DSC) with magnetic 
field induced spin density wave (SDW) order. In the presence of the induced SDW, we find that 
there exists no simple direct relationship between NMR signal rate T^^ and low energy local density 
of states while these two quantities are linearly proportional to each other in a pure DSC. In the 
vortex core region, T^^ on ^^O site may exhibit a double-peak behavior, one sharp and one broad, 
as the temperature is increased to the superconductivity transition temperature Tc, in contrast to a 
single broad peak for a pure DSC. The existence of the sharp peak corresponds to the disappearance 
of the induced SDW above a certain temperature Taf which is assumed to be considerably lower 
than Tc. We also show the differences between T^^ on ^^O and that on ^^Cu as a function of 
lattice site at different temperatures and magnetic fields. Our results obtained from the scenario of 
the vortex with induced SDW is consistent with recent NMR and scanning tunneling microscopy 
experiments. 

PACS numbers: 76.60.-k, 74.60.Ec, 74.25.Ha 



I. INTRODUCTION 



Intensive efforts have been focused recently on the na- 
ture of low-lying excitation spectra around a vortex core 
in high temperature superconductors (HTS). The excita- 
tions around the core play a fundamental role in deter- 
mining physical property of a superconductor. It has now 
been established both experimentall3iiA2i2i^i^ and the- 
oretically^-^-^-^Q-^^-^^-^^-^"^-^^-^^ that a spin density wave 
(SDW) order could be induced and pinned at the vor- 
tex lattice by a strong magnetic field for both under- 
and optimally- doped HTS. In neutron scattering exper- 
iments by Lake et al}^ a remarkable antiferromagnetic 
(AF)-like SDW was observed in optimally- and under- 
doped La2-ajSra;Cu04 under a strong magnetic field. The 
muon spin rotation measurement by Miller et al? stud- 
ied the internal magnetic field distribution in the vortex 
core of underdoped YBa2Cu306+x7 and it was revealed 
a feature in the high-field tail which fits well to a model 
with static alternating magnetization. Recent nuclear 
magnetic resonance (NMR) measurements by Mitrovic 
et al.^ studied the spatially resolved NMR signal in the 
mixed state of optimally doped YBa2Cu307_5, and they 
found strong AF fluctuations outside the cores and rather 
different electronic states inside the vortex cores. An- 
other NMR experiment by Kakuyanagi et al.^ investi- 
gated the magnetism in and around the vortex core of 
nearly optimally-doped Tl2Ba2Cu06+5, the NMR signal 
rate T^^ at Tl site provides a direct evidence that the 
AF spin correlation is significantly enhanced in the vor- 
tex core region. From theoretical point of view, various 
approaches^'^'^'^^'^^ suggest the existence of an induced 
AF order inside the core. Some recent studiesi2ii2iiiiiSii& 
also reveal that the AF order could propagate outside the 



vortex cores and form a SDW. 

As proposed in several articlesiSiiSii^ , that the spa- 
tially imaging NMR experiment can be a powerful tool 
to investigate the exotic electronic structure around the 
vortex cores. The frequency dependence of NMR sig- 
nal rate reflects the internal magnetic field distribution 
and the electronic and magnetic excitations in the mixed 
state can be probed through the temperature and site 
-dependence of T^^ . The spatial variation of the vortex 
lattice in the NMR experiments can be resolved by the 
distribution of internal magnetic field. In the present pa- 
per, we study the NMR theory for HTS with the effect 
of the magnetic field induced SDW being taken into ac- 
count. Our approach is based upon an effective model 
Hamiltonian with a (i-wave pairing interaction Vdsc be- 
tween nearest neighboring sites and a competing onsite 
Coulomb repulsion U which may generate a competing 
AF order. The parameters are chosen in such a way 
that only d-wave superconductivity (DSC) prevails in the 
optimally doped HTS samples when the magnetic field 
B=0. Using the Bogoliubov-de Gennes (BdG) equations. 
The dynamic spin-spin correlation function between site 
i and site j can be numerically obtained. From this cor- 
relation function and the established methods^S, we are 
able to derive the NMR relaxation rates T^^ on ^^O 
and on ^^Cu nuclei, and from which the NMR signals 
as a function of temperature T and site locationS^ can 
be determined. Our calculation is performed for an op- 
timally doped HTS sample with the chosen U so that 
the B induced SDW would appear only below a critical 
temperature Taf which is considerably smaller than the 
superconductivity transition temperature Tc. We find 
that there exists no simple direct relationship between 
the NMR signal rate T^^ and the low energy local density 
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of states (LDOS) while these two quantities are hnearly 
proportional to each other in a pure d-wave supercon- 
ductor {U = 0). In the core region, we find that T^^ ex- 
hibits a double-peak behavior, one sharp and one broad, 
as the temperature is increased to Tc, in contrast to a 
single broad peak for a pure d-wave superconductor. It 
is found that the low temperature peak becomes sharper 
and moves to lower T when Taf becomes lower. This 
result is consistent with the NMR experiments^i^i^. We 
also show the differences between T^^ on ^^O and that 
on ^^Cu as a function of lattice site at different tempera- 
tures and magnetic fields. In general the NMR signal at 
the ^^Cu site is larger than that at ^^O site, and its mag- 
nitude can be quite enhanced at higher magnetic field. In 
section II we will outline our method and the numerical 
scheme for calculating the NMR relaxation rate. In sec- 
tion III, we will present our numerical results and their 
comparison with the experimental measurements. In ad- 
dition we are also going to compare our work with other 
similar theoretically studies iii^fi*^, there the articlesii*^ 
only consider the case for a pure (i-wave superconductor. 
In all these worksiii^fiiSi, the NMR signals are calculated 
for the square lattice sites which may not rigorously re- 
lated to the ^^O and ^^Cu sites as discussed in the present 
paper. In the last section, we will give a summary and 
discussion of our results. 



II. FORMALISM 

Let us begin with a phenomenological model in which 
interactions describing both DSC and SDW order param- 
eters in a two-dimensional lattice are considered. The 
effective Hamiltonian can be written as: 

+ ^(Ay4^4^+/i.c.) , (2.1) 

where c]^ is the electron creation operator, /i is the chem- 
ical potential, and the summation is over the nearest 
neighboring sites. In the presence of magnetic field B 
normal to the two dimensional plane, the hopping in- 
tegral can be expressed as Uj = toexp[i^ J^ ' ^(i*) ' dr] 
for the nearest neighboring sites j), with = /i/2e as 
the superconducting flux quantum. We assume that the 
strength of B is large enough such that it can be regarded 
as uniform. Here we choose the Landau gauge for the vec- 
tor potential A = {—By, 0,0). Since the internal mag- 
netic field induced by the supercurrent around the vortex 
core, as it will be numerically shown later, is so small as 
compared with B, the above assumption should be justi- 
fied. The induced SDW and the DSC orders in our sys- 
tem are respectively defined as Af = U (cJ^qi— cJ^qi) 
and Aij = Voscici^Cji - c/|C^|)/2, where U and Vdsc 
represent the interaction strengths for SDW and DSC. 



The mean- field Hamiltonian (1) can be diagonalized by 
solving the resulting BdG equations self-consistently 

where the single particle Hamiltonian Hij^a = —tij + 
(Urii^ - li)6ij, and n^^ = En Kr = 

tanh ^2^^^, with f{E) as the Fermi distribu- 
tion function and the electron density = + n^|. 
The DSC order parameter at the zth site is Af^ = 
(Af+,^,, + Af_,^,, - - where Ag = 

Aij exp[z^ jJi"*+rj)/2 • dr] and ex^y denotes the unit 
vector along {x,y) direction. The main procedure of 
our self-consistent calculation is summarized as follows. 
For a random set of initial parameters nia- and A^j, 
the Hamiltonian is numerically diagonalized and the ob- 
tained quasiparticle wave functions u^^ and ^ are used 
to calculate the new parameters for the next iteration 
step. The iteration continues until the relative differ- 
ence of order parameter between two consecutive itera- 
tion steps is less than 10~^. 

In the following section, we need to calculate the local 
density of states (LDOS) according to the formula: 

p,{E) = - Y,[\uh?f\E^ -E) + KJ'f'iEr. + E)] , 

n 

(2.3) 

where f{E) = df{E)/dE. The LDOS is proportional to 
the local differential tunneling conductance at low tem- 
perature which can be measured by STM experiments. 
In addition, we also need to calculate the internal mag- 
netic field distribution of the vortex state. The internal 
magnetic field is computed through the Maxwell equa- 
tion, V X Hinti'r) = ^j(r), where the current j(r) is 
calculated as: 

je.„(r,) = 2|e|Im{ti+e.,„iE[<+e.,.</(^") 

n 

+<+e.„<*(l-/(^n))]}- (2-4) 

The nuclear spin relaxation rate T^^ can be obtained 
from the spin-spin correlation function x+,-{h j^^) 
zero energy through the following site dependent func- 

R{iJ) = lmx+,-{iJ,inn^n^iv)/{^/T)\n^o 

xnTf'iE„)5iE„-E„,). (2.5) 

Since all of the experiments ^lii^ are performed at the 
oxygen, copper, or other nuclei in HTS, we need the ex- 
pressions of the NMR relaxation rate at these sites. In 
this paper we shall mainly consider the contribution from 
the planar oxygen (^^O) and copper {^^Cu) nuclei. The 
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oxygen nucleus is located between two Cu nuclei along 
the X or y axis. The internal magnetic field Hint on 
site is averaged over its two nearest neighboring ^^Cu 
sites. The NMR relaxation rate form factor on nuclei ^^O 
and ^^Cu in momentum space is: Fq^x ~ 1 + cos{qx/2) 
and Fcu ^ {cos{qx) + cos{qy))'^. It can be expressed in 
terms of the real-space spin-spin correlations as follows^^: 

l/"TiT = C[i?(i,i) + l/2^i?(i,j)], (2.6) 
1/63TiT = B[R{i,i) + I/2Y, Rihj) 

3 

III 

+l/4^i?(i,i)], (2.7) 

3 

I II III 

here the notation of XI and ^ means summation 
over two nearest-neighboring sites along the x axis, four 
next-neighboring-sites, and four next-next-neighboring 
sites, respectively. Here we calculate T^^ in terms of 
arbitrary unit and choose the ratio of the coefficientsSi 
5/C= 1.583. 

Throughout this paper, we use the following parame- 
ters: U = 2.4, Vdsc = 1-2, the linear dimension of the 
unit cell of the vortex lattice is chosen as Nx x Ny = 
40 X 20 sites and the number of the unit cells Mx x My = 
20 X 40, and the hole doping level x = 0.15. This choice 
corresponds to a uniform magnetic field B 37T. We 
set to = <^ = 1- The standard procedures'^ are followed 
to introduce magnetic unit cells, where each unit cell con- 
tains two superconducting flux quanta. By introducing 
the quasi- momentum of the magnetic Bloch state, we ob- 
tain the wave function under the periodic boundary con- 
dition whose region covers many unit cells. 



III. NUMERICAL RESULTS AND COMPARING 
WITH EXPERIMENTS 

Our previous numerical calculations^' show clearly 
that the induced AF order exists both inside and out- 
side the vortex cores, and behaves like an inhomoge- 
neous SDW with the same wave length in the x and y 
directions. Fig. 1 depicts the SDW amplitude or the 
staggered magnetization M^=Af^^/U distribution in 
a square lattice. It is seen that the induced SDW or- 
der reaches its maximum value at the vortex core center 
(V), zero value at the saddle-point midway between two 
nearest neighboring vortices (S), and in between at the 
center of the squared vortex lattice unit cell (C). Using 
Eq.(2.4), the internal magnetic field Hint at different site 
can be numerically determined. In the upper right inset 
of Fig. 1 we show the histogram of Hint^ where P{Hint) 
measures the number of lattice sites (or plaquettes) with 
fixed Hint- We can approximately identify each site (V, 
S, C) including the sites around it in the vortex lattice 



to the distribution. The correspondence between the site 
position and Hint can roughly be established. It is easy 
to see that the magnitude of Hint is decreasing along the 
path (V^S^C) and its direction is always perpendicular 
to the lattice plane. In a narrow region. Hint may even 
become negative. This is very different from the case for 
a pure DSC where Hint is always positive Our his- 
togram of Hint distribution is consistent with the exper- 
imental dataS except there the region of negative Hint is 
somewhat wider than ours. The experimentally observed 
negative Hint may indicate the existence of magnetic field 
induced SDW around the vortex cores. 

The NMR signal at the maximum cutoff of the his- 
togram as a function of internal magnetic field comes 
from the V-site. The minimum cutoff is from the C-site. 
The peak value of the internal magnetic field comes from 
the S-site. Contrary to the case for a pure d-wave super- 
conductor, the SDW order outside the vortex core may 
have a remarkable contribution to both the internal mag- 
netic distribution and the NMR signal. It is supposed 
that ^^O NMR experiment can provide direct informa- 
tion on the LDOS as the antiferromagnetism spin fluc- 
tuation is cancelled^^'^^. From our study, we find that 
in the presence of SDW vortex cores, the internal mag- 
netic field on ^^O site can not be cancelled exactly and 
appreciable residue exists. 

Next, we study the relationship between low energy 




FIG. 1: The amplitude distribution of the SDW order pa- 
rameter M[ (a) in one magnetic unit cell. V, S, C points 
represents the vortex core center, the saddle point midway 
between two nearest neighboring vortices, and the center of 
the squared vortex lattice unit cell, respectively. The inset 
shows the histogram P{Hint) of the induced magnetic field 
distribution where Hint is in a unit of 10^ Gauss. The calcu- 
lation is performed with = 2.4 and Vdsc = 1-2. 
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FIG. 2: The spatial distribution images of the LDOS at the 
energy ^ = (a) and NMR relaxation rate l/Ti{i)T (b) in a 
vortex unit cell for a pure DSC with [7=0. Other parameter 
values are the same as those in Fig. 1. 



LDOS and NMR signal rate 1/Ti. In a pure d-wave 
superconductor (DSC), a linear relationship is found be- 
tween 1/Ti(z) - R{i,i) and the LDOS pi{E = 0)^^ at 
zero temperature. We have checked our program and 
reproduced similar results. Figure 2 depicts the spatial 
distribution images of LDOS pi{E = 0) (a), and l/Ti(i)T 
(b) for a pure DSC with U =0 at low temperature. In this 
case, the core state shows zero-energy peak^^ and nodal 
directions are along x = y and x = —y. The images 
in Fig. 2(a) and 2(b) are similar to each other and the 
results can be regarded as linearly proportional to each 
other. Thus the NMR signal for a pure DSC can be esti- 
mated quantitatively from its corresponding zero energy 
LDOS according to the linear relationship. 

From Fig. 2(a), the LDOS of quasiparticles decreases 
along the path (V^C^S). In the nodal directions, ap- 
preciable weight of the low-energy states extending from 
the vortex core still remains there. Our numerical re- 
sult indicates that the LDOS of the low-energy states 
at the C-site is slightly larger than that of the S-site, 
although C-site is farther away from the vortex center. 
This is different from the work by Takigawa et al}^ where 
the LDOS at the S site is found to be larger than that 



at the C site. The authors^'' there employed the sym- 
metrical gauge and they found that their unit vector of 
vortex lattice is oriented along 45^ from the Cu-0 bond. 
But our result is consistent with the study of Knapp et 
al. ^^who used a different approach and calculated the 
NMR signal at various sites for a pure DSC. Furthermore, 
a recent small-angle neutron scattering measurement^^! 
for Lai.83Sro.i7Cu04+5 indicates that the unit vector of 
vortex lattice is clearly oriented along the Cu-0 bonds, 
consistent with our result in Fig. 1. The NMR experi- 
ments^'^'^'^ also show that the NMR signal first decreases 
as one moves away from the V-site to S-site and then in- 
creases at the C-site, consistent with the results in Fig. 2. 
Here we notice that the NMR signal at C site is slightly 
larger than that at S-site, their difference could be seen 
in Fig. 6 and Fig. 7. It appears that the above analysis 
could provide a basic understanding of the spatial distri- 
bution of the experimentally observed NMR relaxation 
rates in HTS in the framework of pure DSC. However, 
it is well known that the theoretical predicted LDOS^, 
which shows a strong zero energy peak for a pure DSC at 
the vortex core is in contradictory to the STM measure- 
ments22iS4 on HTS. In order to overcome this difficulty, 
the magnetic induced AF order has been introduced^- 
to account for the experiments. In the following we shall 
investigate the effect of induced AF or SDW order on the 
NMR theory in HTS. 

The presence of the induced SDW order in the mixed 
state may lead to violation of the linear relationship be- 
tween the LDOS at E={) and the NMR relaxation rate 
1/Ti(z) . This can be seen clearly in Fig. 3, where 
we show the spatial images of the LDOS (a) and the 
l/Ti(i)T (b) according to our numerical calculations for 
U=2A. There is little resemblance between Figs. 3(a) 
and 3(b), and the linear relationship does not exists here. 
Fig. 3(a) indicates that the LDOS is smallest at V-site, 
and largest at S-site, it clearly can not be applied to ex- 
plain the NMR experiments^'^'^'^ as in the case for a pure 
DSC . But the LDOS obtained here exhibits no zero en- 
ergy peak at the vortex core which is consistent with the 
STM experiments^^'^^. The averaged LDOS of quasipar- 
ticles in Fig. 3(a) ([/ =2.4) is greatly suppressed than that 
in Fig. 2(a)([/=0). This is because the presence of the 
induced SDW enhances the insulating nature of the sys- 
tem and thus reduces the number of the quasiparticles. 
When [7=2.4, there are two contributions to the NMR 
signal l/Ti(z)T, at the low temperature, the dominant 
one is from the induced SDW order and the minor one is 
from the low energy quasiparticle states. The LDOS is 
the smallest at the vortex core center but the SDW order 
is at its maximum as shown in Fig. 1. From Fig. 3(b), it 
is easy to see that the spatially distributed NMR signal 
l/Ti{i)T which has the largest value at the vortex cen- 
ter because of the induced AF order there and decreases 
along the path (V-^C^S). This feature is consistent with 
what has been observed by the experiments^'^'^'^ and in 
sharp contrast with Takigawa et al^s results^^. They 
claimed that T^^ at V-site is smaller than the neigh- 



5 



(a) LDOS 




5 10 15 20 



FIG. 3: The spatial distribution images of the LDOS at the 
energy ^ = (a) and NMR relaxation rate l/Ti{i)T (b) in a 
vortex unit cell for a DSC with induced SDW. All parameter 
values are the same as those in Fig. 1. 



boring sites and the NMR signal rate at S-site is higher 
than that of C-site. The choice of lattice size in their 
calculations corresponds to a very strong magnetic field 
at about 51 Tesla. Their results seem to be inconsistent 
with the measurement of on YBa2Cu3 07_(5^. 

The NMR signals also exhibit the temperature T and 
magnetic field B dependences. Our approach described 
in Section II is easily extended to study the finite T 
case, but not the B dependence because weaker B im- 
plies larger size calculation which we are not able to do 
at the present moment. Since the observed NMR signals 
are coming from the oxygen, copper, or other atoms in 
HTS, here we shall mainly consider the temperature de- 
pendence of the NMR signal rate from the ^^O 
sites using Eq. (2.6). It needs to be point out that the 
value of R{i^j 7^ i) in Eq. (2.6) is of the same order with 
R{i^i) and its contribution should not be neglected. In 
Fig. 4, the NMR relaxation rates l/Ti = l/^^Ti on ^^O 
site normalized by its value at Tc as a function of 
T are plotted for sites V, C and S. The calculations are 
performed using [7=2.4 except the dashed curve associ- 
ated with the V site where U=0 and is obtained for a 
pure DSC. The inset plots the temperature dependence 
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FIG. 4: Temperature dependences of ^^O NMR signals 1/Ti 
normalized by its value 1/Tf at Tc are plotted as a function of 
T/Tc at the sites V, S, C. The dashed line associated with the 
site V is obtained for a pure DSC by using [7=0. The inset 
plots the temperature dependence of SDW order parameters 
at V-site and DSC order at S-site. 



of the staggered magnetization at V-site and DSC order 
at S-site. Our calculations show that the vast major- 
ity sites outside the vortex core exhibit long relaxation 
time, but few sites around the vortex core exhibit short 
relaxation time, which contains useful information on the 
site-dependent low energy excitations around the vortex. 
The maximum magnitude of ^^T^^ at V-site is about 
two orders of magnitude larger than that far from the 
vortex core, and the NMR signal rate at C-site is also 
larger than that at S-site. Because both NMR signals 
at C and S sites are very weak as compared to that at 
V site, their distinction is not reflected in Fig. 4 below 
T < 0.4Tc. The minimal value of the staggered magneti- 
zation occurs near the saddle point S. It reflects the fact 
that the SDW order at C site is larger than that at S site 
as shown in Fig. 1. This result clearly indicates that the 
presence of the induced SDW order makes remarkable 
contribution to NMR signals. At the vortex core cen- 
ter V, the evolution of T^^ with temperature exhibits a 
double peak structure below Tc; the sharp peak occurs 
at lower T and the broad peak is at higher T. This is in 
contrast to a single peak obtained with U =0 for a pure 
DSC^^, and it can be clearly seen from the dashed curve 
below a critical temperature Taf ^ 0.54Tc and the solid 
curve above Taf associated with the V site in Fig. 4. 
The inset shows that above the temperature Taf, the 
induced staggered magnetization vanishes and only DSC 
order prevails up to Tc. The sharp peak at 0.45Tc is 
below Taf and the broad peak at 0.64Tc is originated 
in the pure DSC order. It is obvious that the value of 
the critical temperature Taf is U or sample dependent, 
and the position of the sharp peak can made to a lower 
T if a smaller U or Taf is chosen for the calculation. In 
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FIG. 5: Site-dependent relaxation rate T^"^ for ^^Cu (a), O 
nuclei of pure d-wave case and ^^Cu (c), O (d) site of SDW 
vortex core case. 



case the value of Taf is slightly higher than the posi- 
tion of the broad peak, only a single peak for the NMR 
signal is predicted. When Taf is close or larger than 
Tc, the NMR signal would be a increasing function of 
T and may not show any peak structure. We also have 
done the same calculation for the NMR relaxation rates 
on ^^Cu nuclei on our model lattice sites (1/Ti(i)), and 
found similar T and site dependence as shown in Fig. 4. 
In the experimental aspect, a sharp peak has been ob- 
served at T - 0.235Tc ^ in TlsBasCuOe while there is 
no observed peak in YBa2Cu307 up to T ~ 0.272Tc ^. 
Therefore the value U in YBa2Cu307 should be larger 
than that in Tl2Ba2Cu06. This conclusion is opposite to 
the theoretical prediction made by Takigawa et al^^ 

We also would like to examine the difference in the spa- 
tial distributions of the NMR signals from ^^Cu and ^^O 
nuclei. In Fig. 5, The site distribution maps of Tf^ for 
^^Cu and for ^^O in a pure DSC {U=0) are respectively 
illustrated in (a) and (b). Both of these maps reach a 
similar peak value at the vortex core and then fall away 
smoothly from it. In the vortex core, the magnitude of 
T^^ at the ^^Cu site is about three times larger than that 
at the ^^O site. In Figs. 5(c) and 5(d), we present the 
results for ^^Cu and ^^O nuclei in a DSC with induced 
SDW ([7=2. 4). The site-dependent spectra for both nu- 
clei show staggered oscillations in the vortex core, there 
the NMR signal rate from the ^^O site is about one fourth 
of the value from that of ^^Cu site. This is because the 
AF order at the ^^O sites could be dramatically weak- 
ened by the partial cancellation of the staggered mag- 
netization between two nearest neighboring sites on the 
original lattice. It is useful to point out that the NMR 
signal ^^T^^ outside the vortex core region is one order 
of magnitude smaller than that of ^^T^~^. In the presence 
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FIG. 6: Site-dependent relaxation rate T^"^ for two different 
temperatures with U—2A and the magnetic field B=37T at 
^^Cu (a) and at ^^O (b). The log plot is shown as inset. 



of the induced SDW, the NMR signal rate is enhanced 
at the core center and suppressed away from the vortex 
core. For both U=0 and [7=2.4 cases, T^^ reaches a 
maximum at the core center while its value at the C-site 
is always larger than that at the S-site. 

Next, the temperature dependence of NMR signal rates 
for the case of U =2.4 is examined. Here, we plot the vari- 
ation of T^^ along the path V ^ S ^ C for ^^Cu nuclei 
in Fig. 6(a) and for ^^O nuclei Fig. 6(b) at two different 
temperatures. The values of T^^ in Figs. 6(a) and 6(b) 
exhibit similar behavior and oscillate from site to site. 
The amplitude of the oscillation reflects the underlying 
staggered magnetization and is apparently weakened at 
higher temperature. In fact, the T dependence of NMR 
signal rate at three typical points has been seen from Fig. 
3. Figure 6 shows more clearly that NMR signal rate 
at each site increases with temperature and its value at 
site C is always larger than that at site S. We have re- 
produced the power law relation T^^ ~ for a pure 
DSC^^ at zero magnetic field and at low temperature. 



7 




V s c 



0J2 




FIG. 7: Site-dependent relaxation time T^^ for two different 
magnetic fields with [7=2.4 and the temperature T/Tc=0.35 
at ^^Cu (a) and at ^^O (b). The log plot is shown as inset. 



In the presence of vortices with induced SDW order, it 
should not be surprised that T^^ does not follow this 
law. 

Finally, we show the site dependence of NMR signal 
rates on ^^Cu and on ^^O nuclei at two different mag- 
netic field B respectively in Figs. 7(a) and 7(b). The 
NMR signals T^^ in Figs. 7(a) and 7(b) are similar to 
each other. Since the magnitude of the induced SDW 
order is strongly dependent on the external magnetic 
field B — higher B leads to more pronounced SDW or- 
der which may in turn yield a larger NMR signal. This 
conclusion is consistent with the NMR experiments^. As 
a result, the amplitude of oscillation at different sites, 
which measures the strength of the induced staggered 
magnetization, is dramatically reduced when the mag- 
netic field B is decreased from 37T to 27T as shown in 
Fig. 7. Here B = 27T corresponds to a vortex lattice 
unit cell of 48 x 24 sites. On the other hand, the NMR 
signal rate at V-site is suppressed with increasing B for 
a pure DSC — which reflects the density of quasiparticle 



around the vortex core is reduced due to the enhance- 
ment of quasiparticle leakage along the nodal directions. 
It was argued by Takigawa et ali^i that when AF order 
inside the vortex core is large enough, T^^ becomes the 
absolute minimum at V-site. This conjecture^^ is based 
on the assumption that the linear relation between T^^ 
and LDOS still exists even in the presence of the SDW 
order. On the other hand, our results indicate strongly 
that with the increasing B, the enhancement of AF order 
inside the vortex core may lead to more pronounced T^^ 
value at the V-site. This is very different from what has 
been predicted in Ref. 31. Finally we point out that ex- 
isting experimental data^« show the presence of apprecia- 
ble NMR signal rate at S-site and its value is insensitive 
to T. But in our calculations there is a very weak NMR 
signal at S-site. The discrepancy between experimental 
data and our results may be due to the reason unclear to 
us at this moment. 



IV. SUMMARY 

Based on a model Hamiltonian as described in Sections 
I and II, we have investigated the site- and temperature- 
dependence of NMR signal rate T^^ near the vortex core 
in an optimally doped HTS. The calculations are carried 
out with {U = 2.4) and without {U = 0) the magnetic 
field induced SDW order. Although the relative values 
of T^^ at V, S, and C sites near the vortex core ob- 
tained theoretically for a pure DSC are consistent with 
recent experiments ^^^i^i^ , the temperature dependent of 
T^^ thus derived fails to explain the sharp peak observed 
in NMR experiments^ for Tl2Ba2Cu06+5 at low temper- 
ature. Moreover, the LDOS obtained for a pure DSC 
— has a zero energy peak at the core center and this is 
in contradictory to the STM observations^^' All these 
difficulties could somewhat be removed by introducing 
the field induced SDW pinned at the vortex cores as it 
has been studied in the present paper. With the induced 
SDW order, the LDOS exhibits no zero energy peak at 
the core center and is consistent with experiments. Un- 
like the case for a pure DSC, there exists no linear rela- 
tionship between zero energy LDOS and T^^. We show 
that SDW order could strongly enhance the NMR sig- 
nal rate, and as a result, the NMR signal T^^ from ^^O 
site has its largest value at the core center V, and small- 
est value at the saddle point S site (see Figs. 6 and 7). 
The temperature dependent of T^^ at V site exhibits 
a sharp peak at low temperature which is in agreement 
with the NMR experiments^'^. Since the strength of the 
SDW order depends on the magnitude of B, the value of 
T^^ should also be enhanced by increasing B. We also 
compare T^^ from ^^Cu and ^^O nuclei at different T 
and B, their essential features are similar to each other. 
The magnitude of ^^T^^ outside the vortex core region is 
about one order of magnitude smaller than that of ^^T^^ 
while these two values are closer to each other inside the 
vortex core. Finally, we would like to point out that 
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the difference between our results and those of a similar 
work^^ has been discussed in detail in Section III. An- 
other fundamental difference is that their LDOS^^ at the 
vortex center for a pure DSC {U = 0) shows double peaks 
near E=0 while only a single peak is obtained in previ- 
ous^^and present works. It is our hope that the present 
calculation may shed more light on the understanding of 
recent NMR and STM experiments for the mixed state of 
HTS in a unified picture, and stimulate more theoretical 
activity in this field. 
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